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Summary. — A computer code to simulate temporal evolution of overdensities 
in a one- dimensional Universe is presented for didactic purposes. The formation of 
large scale structures in this one-dimensional universe can be studied both in matter 
or radiation dominated eras. Since large scale structures are already observed at 
z>7, primordial dark matter overdensities 5dm which are 90 times larger than the 
observed barionic 5_g in the cosmic microwave background are required at z ~ 1000. 
This makes possible non-linear gravitational collapse at redshift z > 7 and the 
formation of the structures. Primordial perturbations 5b ~ 10~^ do not leave the 
linear regime of growth without the aid of dark matter's potential wells. This code 
is suitable for commercial worksheets like MSExcel, StarOfEce, or OpenOfRce. 

PACS 98.80.H - Cosmology. 
PACS 95.35. -hd - Dark Matter. 
PACS 98.65.-r - Large Scale Structure. 



1. Introduction 

Computer codes for simulating the evolution of large scale structure have been mas- 
sively used in cosmology since years 60s There exist complex codes to study nu- 
merically the evolution of density perturbations in the expanding Universe E] . In 
those codes it is possible to tune several cosmological parameters representing different 
physical conditions like coupling between barions and photons, dependence of density pa- 
rameter VLq on the abundance of various components (hot and cold dark matter, barions, 
cosmological constant). 

The code here presented is based on the Friedmann solution (1922) for the Euclidean 
Universe as expanding background. Moreover the evolution of a density perturbation 
does not affect the expansion itself. 

The limitation to a one-dimensional analysis greatly reduces mathematical complexity 
while allows to obtain qualitative and quantitative predictions on the density perturba- 
tion evolution. 

The evolution of density perturbations of amplitude 5 = 5p/ p = 10^^ starting from 
the decoupling time down to present time is calculated. The requirement of dark matter 
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existence in order to explain the formation of large scale structures at redshift z > 7 is 
shown. 

The advantage of one-dimensional model is the easiness of computation even with a 
worksheet like MSExcel, it is indicated in teaching to high school up to undergraduate 
students. The calculation of equations' coefficients is the only remaining difficulty. 

2. Overdensities in the Universe: a brief history 

Friedmann solution of Einstein equations in the matter era (whose state equation is 
P = 0) yields a{t) ~ t^/a^ 

This regime occurs after the so-called equidensity time in which p of matter equals p 
of radiation at z = 40.000 i.e. 7 years after the Big Bang. 

Going backwards in time, between matter era and radiation era, which is reproducible 
in the computer code changing the expansion rate according to a{t) t^/^, there is 
Meszaros effect: density perturbations do not grow because free fall time is larger than 
expansion time, they stay frozen. 

In particular the derivative of a{t) goes to infinity near singularity, therefore in the 
early stages of Universe's history, gravitational collapse cannot occur because Universe 
expands more quickly than any other dynamical time. 

3. Matter era 

In the matter era we can study the case of perturbations larger than cosmological 
horizon (the Hubble radius) and the cases of linear growth and non-linear collapse. 

In this one-dimensional model, overdensities S are represented as two pointlike masses 
m separated by a distance d such that A/SwcPp = m, where rho is the average density 
of the Universe at decoupling time z = 1000. The volume of the sphere of radius d 
contains the mass m -considered pointlike- at the average density of the Universe at the 
decoupling time. 

In order to simulate an overdensity of (5 = 10~^ it is sufficient to multiply one of the 
masses to by (1 -I- 10~^). 

4. Linear and nonUnear growth of overdensities 

Whatever is the size of the overdensity, it grows according to the approximate law 
log{5{t)) ~ log{a{t)) =^ S{t) oc a{t) 

We show this behaviour from globular cluster size (IO^Mq) to superclusters (IO^^Mq), 
in figure 1. 

Nonlinear growth before z « 7 is necessary for the onset of large scale structures. 
In figure 2 dark matter overdensities 90 times larger than the observed at decoupling 
time in cosmic background radiation (CBR) are capable to collapse before present time, 
regardless on their initial mass or size. 

5. Computer code 

In order to reduce the number if iteractions, and at the same time, maintain enough 
temporal resolution at early stages of Universe's history, we adopted an exponential 
stepping, ti oc e*-'' '-'. This exponential time is used for computing the discrete expansion 
factor d{ti . . . ijv) according to the opportune dynamical regime. Once ready the vector 
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Fig. 1. - Linear growth of density perturbations with 5o = 10~^ equal to those observed in CBR. 




Fig. 2. - Non-linear growth of perturbations: density perturbations evolution of So = 90 times 
those observed in CBR up to their non linear collapse. 
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Gravitation in one- dim expanding Universe 




Fig. 3. - Scheme of one-dimensional Universe's computer code 



a(t), it is used in the calculation of the gravitational force between the two masses, at 
the end of each iteraction, in order to simulate the cosmological expansion. This code 
has been implemented in a MSExcel worksheet. 

* * * 

This work has been developped at Technical Industrial Institute "Giuseppe Armellini" 
of Rome during tlu^ course^ "Lcggcrc il CicJo" -2003 edition- held by the author under the 
auspices of the Italian Ministere of Education. 
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